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The spin of an electron trapped in a quantum dot is a promising candidate implementation of a 
qubit for quantum information processing. We study the central spin problem of the effect of the 
hyperfine interaction between such an electron and a large number of nuclear moments. Using a 
spin coherent path integral, we show that in this limit the electron spin evolution is well described 
by classical dynamics of both the nuclear and electron spins. We then introduce approximate yet 
systematic methods to analyze aspects of the classical dynamics, and discuss the importance of the 
exact integrability of the central spin Hamiltonian. This is compared with numerical simulation. 
Finally, we obtain the asymptotic long time decay of the electron spin polarization. We show that 
this is insensitive to integrability, and determined instead by the transfer of angular momentum to 
very weakly coupled spins far from the center of the quantum dot. The specific form of the decay 
is shown to depend sensitively on the form of the electronic wavefunction. 

PACS numbers: 73.21. La, 85. 35. Be, 72.25.Rb 



I. INTRODUCTION 

A single electron in a quantum dot provides a ver- 
satile experimental realization of a single quantum bit 
("qubit"). The electron can be trapped with a high de- 
gree of certainty by utilizing the Coulomb blockade effect, 
in which a large charging energy prevents the motion of 
electrons on or off the dot. The two level degree of free- 
dom of the qubit is provided then by the electron spin. 
The electron can be readily manipulated electrically and 
optically. A limiting factor in such manipulations is "de- 
coherence" , caused by coupling of the electron spin to 
other degrees of freedom, fn high-quality GaAs quantum 
dots, in which all extrinsic sources of decoherence can be 
removed, the hyperfine interaction with nuclear moments 
becomes dominant 

In a typical quantum dot, the electron spin is coupled 
to a very large number N ~ 10 4 — 10 6 of nuclear moments, 
with a range of hyperfine exchange constants determined 
by the wavefunction of the confined electron. This varies 
smoothly from near zero at the fringes of the dot to a 
maximum value near the dot center. To a rather good 
approximation, there is no direct interaction between dif- 
ferent nuclear moments. The latter is present only due to 
dipolar coupling, but is very small due to the tiny nuclear 
dipole moment. When it is neglected (which is expected 
to be a good approximation for a relatively long time 
~ 10~ 3 sec), the Hamiltonian Tt of the electron-nuclear 
system has the "central spin" form, in which each nuclear 
spin operator Ii (of spin /) is coupled only to the central 
electron spin (1/2) operator S: 

H = J2 a ^ ■ $ ■ (!) 

i 

The exchange constants a* vary with the amplitude of 
the electron wavefunction in the Quantum dot, at the 
positions ?i of the nuclei: a, ~ | -0 (7^ ) | 2 . 



The central spin model is therefore of considerable ex- 
perimental relevance. It is also of particular theoreti- 
cal interest. Most intriguingly, it is known that elec- 
tron spin decoherence within this model is highly non- 
conventional. For instance, the relaxation of an initially 
polarized electron spin is non-exponential. Indeed, nu- 
merical simulations for some specific forms of the electron 
wavefunction indicated an inverse-logarithmic remanent 
polarization at long times^^ 

with the initial condition S z (0) = +1/2 and random un- 
polarized initial conditions for the nuclear spin state (see 
Sec. Ill Ap ; these numerical results suggested a w 1 . 

One possible reason to expect unusual behavior in the 
central spin problem is that it is integrable. In fact, the 
central spin Hamiltonian is a member of a large family of 
models studied by Gaudin, possessing an infinite number 
of commuting constants of the motion, both classically 
and quantum mechanically. The Gaudin models unify 
the central spin problem with another problem of recent 
interest, the reduced BCS Hamiltonian for a small disor- 
dered superconducting grain. This integrability has at- 
tracted additional theoretical attention;^ but as is often 
the case, it is quite difficult to extract real-time dynam- 
ical properties despite the model's formal solubility. In 
any case, it is not clear from these studies whether the 
unconventional long-time dynamics should be ascribed to 
integrability or has another origin. 

Another controversial aspect of the central spin prob- 
lem is the significance of quantum coherence to the ob- 
served behavior. The pioneering studies in Refs. Si! 
insist upon a full quantum treatment, but employ ap- 
proximations valid only in a large applied magnetic field 
or for nearly full nuclear polarization. The work of Reff^ 
instead argued on physical grounds for the validity of 
the classical approximation, and integrated the classical 
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equations of motion for the spins numerically, with how- 
ever no assumptions on the magnetic field or polarization 
(the focus was largely on the low field and polarization 
situation most relevant to experiment). Refsl4U5lll0l argue 
for a hybrid approach in which effective classical equa- 
tions of motion are employed using renormalized param- 
eters, finding good agreement with fully quantum simu- 
lations on small systems of 15 — 20 spins. 

In this paper, we will attempt to clarify the above is- 
sues through a combination of techniques. We will first 
argue, in agreement with Ref,6, that the classical approx- 
imation for the nuclear spins is indeed justified for mea- 
surements of the electron spin dynamics in the physically 
relevant large N limit. We show this by a simple direct 
analysis of the exact quantum expression ( "Keldysh-like 
form") for the time-dependent expectation value of the 
electron spin, thereby obtaining the explicit connection 
between the properly defined quantum observable and 
the classical dynamics. The regime of validity of this 
approximation is discussed, and open issues highlighted. 

Next, having established the correctness of the classical 
limit, we attack the problem of the anomalous long-time 
spin dynamics. Despite the simplification of the clas- 
sical limit, the dynamics remains highly non-linear and 
non-trivial, and still is not easily accessed using integra- 
bility - though some progress has been made in this way 
for the BCS probleml&ii. We will argue that the ob- 
served long-time behavior is due not to integrability but 
rather to the effects of very weak but non-zero coupling 
of the electron to spins far from the center of the dot, in 
the tails of the electron wavefunction. As time proceeds, 
successively more weakly coupled nuclear spins are able 
to influence the electron spin dynamics, leading to a slow 
change in the latter's behavior. To address the effects of 
integrability, we employ a phase space approach inspired 
by fundamentals of statistical mechanics. In general, one 
expects (the ergodicity assumption) a Hamiltonian sys- 
tem to equally explore over sufficiently long times all of 
its phase space accessible by the constraints of its conser- 
vation laws. An important feature particular to the cen- 
tral spin problem is the presence of very weakly coupled 
spins in the tails of the electron wave function. Even for 
very long times, some such spins have not yet been able 
to evolve appreciably, further restricting the phase space 
which has been accessed up to this time. We capture 
this physics through a crude but effective approximation. 
For any given time t, we divide the nuclear spins into two 
groups: those with sufficiently large couplings to the elec- 
tron spin to allow them to precess through an angle of 
0(1), and those with weaker couplings that have not. We 
then neglect completely the latter "decoupled" spins, and 
treat the former "coupled" spins as an isolated subsys- 
tem, and assume that it has fully explored its accessible 
phase space already in the time t. In this "instantaneous 
ergodicity" approximation, the long time dynamics is due 
to the continual growth of the ergodic subsystem by in- 
clusion of additional weakly-coupled spins as time pro- 
gresses. 



To study the properties of the ergodic average of the 
coupled subsystem at any given time, one must still deal 
with the effects of integrability. Specifically, since the 
central spin problem is integrable for any number of nu- 
clear spins, the accessible phase space is especially highly 
constrained due to the infinite number of constants of the 
motion. We contend, however, that the constraints of in- 
tegrability are secondary to the already-described physics 
of weakly coupled spins. The basis of this conclusion is 
a re-organization of the integrals of motion into an or- 
dered set in which successive integrals manifestly have 
less influence upon the electron spin. The first four (and 
therefore most important) integrals of motion in this set 
are the three components of the total angular momen- 
tum and the energy. We find that ignoring all remaining 
integrals of motion, i.e. replacing averages over times 
less than t by a uniform phase space average in the cou- 
pled system at time t without these constraints, already 
reproduces the observed non-exponential long-time spin 
relaxation. We show how additional integrals of motion 
can be systematically included in this analysis, and com- 
pare the results to numerical simulations. 

For the specific forms of wavefunction employed in 
prior numerical simulations, an inverse logarithmic relax- 
ation is obtained for long times, in agreement with those 
numerical results. This form is, however, not univer- 
sal, and depends in particular upon the manner in which 
the electron's wavefunction approaches zero. The result 
can be understood rather simply as follows. Because the 
electron is coupled to many nuclei, its characteristic pre- 
cession frequency is much larger than the rate at which 
any of the individual nuclei themselves evolve, as they 
are coupled only to the single electron spin. Under this 
condition, as pointed out in Ref. @, the time-averaged 
electron spin simply follows the evolution of the net hy- 
perfine field 

Hv »' (3) 

j 

Because the hyperfine couplings dj decay with the dis- 
tance of the spin Ij from the center of the dot, the above 
hyperfine field can be crudely approximated by including 
only those N s spins whose couplings are approximately 
constant and equal to the maximal a max : 

N s 

H N Q max y^Ij. (4) 

j=i 

Thus the hyperfine field is approximately proportional 
to the total angular momentum of these N s strongest- 
coupled spins. The total angular momentum of all spins 
is however conserved, and therefore cannot decay but can 
only be redistributed. Because there is no preference for 
how this angular momentum is spread amongst the sys- 
tem's spins, it is natural to expect that over long times, 
the total angular momentum of all N(t) spins which are 
appreciably involved in the dynamics is uniformly dis- 
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tributed amongst them, and we have simply 



N(ty 



(5) 



The time-dependence of N(t) is determined from the fact 
that, at a large time t, the dynamical spins are those 
which have had time to precess, i.e. for which ajSt > 1. 

We can in this way readily understand the logarith- 
mic behavior of Eq. ([2]) for the wavefunctions studied in 
RefsHH. For these Gaussian and exponential wavefunc- 
tions, 



(6) 



where R s gives the "radius" of the dot and determines 
the number of strongly coupled nuclear spins, N s ~ Rf, 
for a d-dimensional dot. The parameter 7 = 1,2 for ex- 
ponential and Gaussian wavefunctions, respectively. The 
dynamical spins at time t are then those inside some ra- 
dius R(t) such that 



R, 



(lni) 



1/7 



(7) 



From Eq.©, we then obtain Eq.© with a = d/j. For 
the two-dimensional quantum dots (d = 2) studied in 
Refsii, this gives a = 2, 1 for the exponential and 
Gaussian wavefunctions, respectively. In Sec|V]we obtain 
these and some other results from the above-described 
analytic calculations. 

The remaining parts of the paper are organized as fol- 
lows. In Sec. UH we use spin coherent path integral to 
derive the semi-classical equation of motion for both elec- 
tron spin and nuclear spins in the large TV (number of 
nuclear spins) limit. In Sec. IIII1 we make a short-time 
approximation to the classical equations of motion and 
then focus on the nuclear subsystem. Having system- 
atically argued the relative importance of the integrals 
of motion, we derive the long time average of effective 
nuclear field by only considering the most important in- 
tegrals of motion. The results from the numerical simu- 
lation and approximate predictions given in Sec. IIIII are 
compared in Sec. lIVI Finally in Sec.fV] with the approxi- 
mation explained above, the asymptotic time dependence 
of (S z (t)) is obtained for gaussian, exponentially decay- 
ing and hard wall confining boundary coupling profiles. 
We conclude in Sec lVIl with a summary and discussion 
of future prospects for theoretical development and com- 
parison with experiment. 



II. MODEL AND SEMICLASSICAL LIMIT 

A. Formulation 

We begin by formulating a specific dynamical prob- 
lem described by the central spin Hamiltonian in Eq. (fTJ) . 



Specifically, we consider a general initial state in which 
the electron spin is polarized up: 

l*(0))= Y. CW.-.mJmi,--- ,m N ) I \^) s , (8) 



where subindex I,S mean nuclear subsystem and elec- 
tron, respectively. And J? = m,. The complex coeffi- 
cients C mij ... >mN will be taken as independent zero mean 
Gaussian random variables for each set of rrii . The dis- 
tribution is then fully specified by the variance 



mi,- ,m N ^m. 



n ' N = P{mi , • • • , m w )S TOiOT / • • • 5 mNm , N , 

(9) 

where p(mi,-- - , Jnjv) is a normalized "classical" prob- 
ability distribution for the z-components of the nuclear 
spins. We will work with the general form 



p(mi,- 



,m N ) = -e h ^ 



(10) 



where Z = (sinh[(I + l/2)h}/ sinh(h/2)) N by normal- 
ization. The parameter h specifies the degree of initial 
nuclear polarization along the +z axis, with h = being 
unpolarized and h — > oo for full polarization. 

The remanent electron spin polarization is then defined 

as 



C(t)^(S*(t)) = (y(0)\S*(t)\y(0)), 
where as usual, in the Heisenberg picture, 



S(t) 



jnt 



Se 



-mt 



(11) 



(12) 



It has been observed numerically 5 that (S(t)) defined 
in this way is essentially self-averaging, i.e independent 
of the realization of the coefficients in Eq.©. In Ap- 
pendix. [3] we demonstrate that this is indeed the case 
for N ^> 1, by showing that the normalized fluctuations 
become exponentially small in this limit: 



(SW-((S'(*))) 2 / tanh(lfe) 



(«*(*)> 



Vtanh[(J+ \)h\ 



N 



(13) 



since the quantity in the brackets is always less than 
unity. Therefore one may approximate 



C(i)«C(() = (S'(t)). 



(14) 



The latter averaged quantity is considerably simpler to 
treat theoretically. 

Using Eqs. (|8I9I10|) one obtains therefore 

C(t) w i (|| Tri[e-' ii ' a e™ t S'e- tHt e- ii ' a ] ||> s , (15) 



with the Boltzmann-like factor 



(16) 
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The partial trace and its contents in Eg. ([15]) are ex- 
tremely reminiscent of a conventional equilibrium cor- 
relation function. The only difference is the presence of 
the electron spin operator inside the Hamiltonian 7i. 

To bring out this connection, we introduce a partial 
Trotter decomposition for the exponentials. We make 
use of the spin coherent states, with the general form for 
spin J: 



\h)=n{e,ct>)\j,+j). 



(17) 



Here the 71(9, <fi) a the rotation operator that changes the 
quantization axis from z to the direction defined by the 
angles (6, 4>). The resolution identity in this basis is 



2J+1 
An 



dh\h) (h\ = 1, 



(18) 



where the integration is over the directions on the unit 
sphere. Similarly the trace is given by 



TrO = 



2J+ 1 

47T 



dh(h\0\h). 



(19) 



We first use these coherent states for the electron spin 
only to decompose the exponentials of ±iTtt, and obtain 
the path-integral expressions 

„-mt 

[dh+ {t)]e lSB [™+lf +e~ jS So dt ' A + (*')•« \h+ (i)) (h+ (0) I , 
[dh-{t)]e- tSB[fl - 

h = ajij, 



iS S° dt,fi -^- a \h-(0))(n-(t)\, 

(20) 



where T t + ( T 4 _ ) is the time-ordering (anti-time-ordering) 
operator. Here and in the following, we will use square 
brackets to denote a functional integration measure. The 
overlap between coherent states in the Trotterization in- 
duces as usual a Berry phase, 



S B [n] 



S I dt' (j) cost 
'o 



(21) 



where h = (sin # cos </>, sin sin 0, cos 0) and S = 1/2. Be- 
cause the different nuclear spin operators in each term 
of the sum defining H (Eq. ([20|) ) commute, the opera- 
tor exponentials can be factored. Using then Eq. (f2"0"]) in 
Eq. (|15[) . and expressing the trace over the nuclear sub- 
spaces using Eq. (fT9"]) . one obtains 



C(t) = J PJcMi J [dh+dh-] 

e i{SB\M-S B [*A) (fi_(t)\S*\h + (t)) 

n<^|e^ /2 C/-(0,t)(7+(0,t)e 



(22) 



hI °' 2 \h„ 



>. 



where Z = Z(Aw/(2I + 1)) N , \hi) is a coherent state for 
the nuclear spin I iy and 



(23) 



Inspecting Eq. (p2]) . we note that the integrand of the 
functional integrals over h±(t) contains a product of N 
factors - expectation values in the state \hi) - one for 
each individual nuclear subspace. Each such factor is a 
functional of h±(t), and is distinguished from the others 
only by its value of a, (and the state rlj)). Since there 
are many (O(N)) such factors all with a* of the same 
order, this product is sharply-peaked/rapidly-oscillating 
and renders the path integral suitable for a saddle point 
expansion (formally it can be rewritten in the usual 
form as the exponential of a sum of N smoothly- varying 
terms). We will see below that the Berry phase terms 
are also rapidly oscillating. 

The saddle point equations are: 



8h±(t') 
± iS B [n± 



E 



ln<f» i |e w f/ 2 Ur(0,t)U < + (0,t)e w "/ 2 |n i > 



= . 



(24) 



This saddle point equation actually leads to the classical 
equations of motion. To see this, we evaluate the func- 
tional derivatives in Eq. (|2~4f . The variation of the Berry 
phase term is well knowni 2 -: 



5h 



= S h± x n±. 



± 



(25) 



The variation of the expectation value is more involved: 

^^ln<n i | e w ?/ 2 C7r(o,t)C/+(0,t)e w ?/2|^. ) 
on±(t') 

= ^iaiS(Ti)f t , , 



{h)t, = 



(n i |e^/ 2 ^-(0,Qt/+(f,t)/>+(0,t')e"f/ 2 |n I ) 

(hi\e hI t/ 2 U~(0,t)U+(0,t)e hI ?/ 2 \hi) 
(fn\e hI !/ 2 Ur(Q, t)&r(if, t)U+(0, t')e hI ^ 2 \h t ) 



(hi\e hI ^ 2 Ur(Q,t) U+ (0, t)e w f I 2 \ hi 



(26) 



Thus the saddle point equations become: 

h±(t') x h±{t') = 53o,(^)± . (27) 

i 

These equations are solved by the restricted condi- 
tion h + (t') = h-(t') = h(t'). We then accordingly 
write U^(t\,t2) = Ui(ti,t 2 ) and one has U~(t 1 ,t 2 ) = 
[Ui(ti, £2)] ■ One obtains 



(28) 





e hI i/ 2 [U i (0,t')]^l i U i (0,i/)e hI '/ 2 


hi) 


(hi 


e hI t 


hi) 
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The "classical" nuclear spin variable denned in this way 
has the form of the expectation value of a standard 
Heisenberg operator. Thus it obeys 



d t li = aiS x Ij, 



with 



S = Sh, 



and the initial condition 
if'»(t = 0) = / 



If(t = 0) = I 



cosh (/i/2) + sinh {h/2)n z i 
cosh (h/2)n\ + sinh (h/2) 
cosh(/i/2) + sinh(/i/2)nf ' 



(29) 



(30) 



(31) 



and If (0) = I 2 , in the unpolarized case when h = 0, the 
above expression reduces to a simple form 1^(0) = I hi. 
Combining the results of Eqs. ([27|) .([28 | and the definition 
in Eq. (|30|) . the saddle point equation itself reduces to 



d t S = H„xS, 



(32) 
(33) 



Because of the polarized initial electron spin state, one 
has 



S(t = 0) = Sz. 



(34) 



Note that Eqs. (|29l31l32l34p are precisely the classical 
equations of motion for this system, with the specified 
initial conditions. The remanent polarization is obtained 
from these classical variables by evaluating the path in- 
tegrals over h±(t') at the saddle point, and performing 
the remaining integrations over the hi variables: 



C(t) = [S'(t;{ni})]. 



(35) 



Here S z (t; {hi}) is the solution to the classical equations 
of motion for the nuclear initial conditions specified in 
Eq. ([3"Tj) . evaluated at time t, and the square brackets 
indicate an average over the nuclear initial conditions: 



[A({ni})} = 



/ I], dhi[cosh (h/2) + sinh (h/2)n z } 21 A({hi}) 



J Hi dhi[cosh (h/2) + sinh (h/2)n 



zl27 



The normalization is determined by the requirement that 
(7(0) = 1/2. 



III. CLASSICAL ANALYSIS 

A. Reduced nuclear equations 

In the previous section we showed the classical equa- 
tions of motion are a good approximation for the dynam- 
ics of the system. In this section we proceed to analyze 
them. The classical equations of motion read 



8 t S = H N x S 
dtlj = OjS x Ij 



(37) 



where Hjv = j a j^-j ac t s as an effective magnetic field, 
which the electron spin precesses about. Here |S| = S = 
1/2 and |Lj| = /. 

Since the magnitude of the effective field Hat = |H/v| 
is of order ~ v~N, while an individual nuclear spin cou- 
ples with only the single electron spin, the timescale of 
the dynamics for the nuclear spins is much longer than 
the electron spin^ As a result, electron spin precesses 
rapidly about an axis which itself moves only much more 
slowly with time: 



S(t) = (S) ( + VS 2 -\(S) t \ 2 Re [( ei + ib^ 1 } , (38) 

where ei , &2 are orthogonal unit vectors in the plane per- 
pendicular to the short-time average of the electron spin 
<S>*, 



S cos 0Qih, 



(39) 



which itself is along the axis rh of the effective nuclear 
field 



m 



Htv 



(40) 



The vectors m, ei,e2 form an orthonormal frame, and 
vary on timescales much slower than the electron pre- 
cession quasi-period 1/Hjv- The parameter #o is the 
initial angle between the electron spin polarization and 
Hjy(t = 0). To study the dynamics of the nuclear spin, it 
is sufficient (for large N) to replace the electron spin by 
its time average, leading to reduced equations of motion 

dtlj = gjth x Ij, (41) 

With gj — S cos 6q ctj . Eqs. (|4"Tj) eliminate the electron 
spin completely, and form a closed description of the nu- 
clear dynamics. They are equivalent to the Hamiltonian 

Hcs = Scos(6 )H N , (42) 

which involves only the nuclear spins. 

B. Conserved quantities 

From the Hamiltonian in Eq. (|4"Tj) . it is obvious the 
magnitude of the effective field Hat is a conserved quan- 
tity. By spin-rotational invariance, the total spin angu- 
lar momentum is conserved, which in the above time- 
averaged electron approximation becomes conservation 
of just the total nuclear spin I = J^i Ij (t ne electron's 
contribution to the total spin of the system is negligi- 
ble for large N). It is readily verified that I is indeed 
conserved directly from Eq. (|4"Tj) . 

In fact one can find a series of conserved quantities for 
the reduced equations of motion, obtained in the same 
way from those for the complete central spin problem 11 



^ 9 = E 



2 9]9q 
9j ~ 9q 



(43) 
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where q = 1, . . . , N runs over all the nuclear spins, and 
we have assumed that no two nuclear spins have the same 
coupling gj to the electron spin (this is a simplifying but 
inessential assumption for our purposes). 

It will be convenient to consider a different linear com- 
bination of these constants of the motion, which are nat- 
urally also conserved: 

C n = Y,9«n q . (44) 
q 

The C n can be defined for any n, and any set containing N 
members of C n 's is completely equivalent to the original 
set of {H q }, because the matrix determinant of this linear 
tranformation from {Ti. q } to {C n } is nonvanishing, as all 
the couplings are different. Also notice that the original 
set of Tl q is linearly independent except for one constraint 

n 

c = J2 n « = °> ( 45 ) 

9=1 

which follows from Eq. (|43"|) . The energy S cos (8o)Hn and 
magnitude of the total spin are contained in this set via 

d = -S 2 cos 2 (8 Q )H 2 N +I 2 Y,9l (46) 

3 

C_i = I 2 - I 2 N, (47) 
where I = |l|. 

C. Collective variables and steady-state relations 

It is useful to define a set of vector quantities, which 
serve as collective variables: 

i', ( 48 ) 

3 

And hi = S cos (9q)Hn, ho = I. The values of h„ deter- 
mine the constants of motion C n (but not vice- versa) : 

n 

n>0 , C n = - £ h fc • h n+i-fc + nl2 ( 49 ) 
fe=i 

—n—X 

n < , C n = h * ■ h »+i+t + nj2 Yl 5" +1 ( 5 °) 

fe=0 j 

In addition, the h„'s obey simple dynamical equations 
obtained by taking moments of the equations of motion, 
Eq.gU): 

d t h n = m x h n+ i . (51) 

Here, we have expressed the spin dynamics completely 
in terms of these collective variables. Any set of N con- 
tiguous collective variables, {hj,, • • • , hfc + jv-i} for & n y k, 
is equivalent to the original spin variables. This is be- 
cause Eq. (l4"5|) defines a linear transformation of the van 



der Monde form, whose determinant is therefore non- 
vanishing provided none of the gj are equal, as assumed. 
Furthermore, the field hk+N appearing in the equation of 
motion (Eq. (|51[) ) for dthk+N-i can be eliminated in favor 
of this set since it is not linearly independent of them, 
leading to closed equations of motion for the collective 
variables. 

The point of view in which the h„ are taken as our 
dynamical variables is useful for elucidating the relative 
importance of the C n for positive n. The constants of 
motion of course prevent the collective variables from 
crgodically exploring all of phase space. We are inter- 
ested primarily in the nuclear hyperfine field hi. From 
Eq. , clearly all C n for n > involve hi , and thereby 
restrict somewhat its evolution. However, the invariants 
for increasing n involve progressively more terms which 
do not include hi. Thus we expect the corresponding 
constraints to be less restrictive on the average nuclear 
field. 

Conversely, to understand the (in)significance of in- 
variants C n with negative n is more readily seen in the 
original nuclear spin variables. Specifically, when n < 0, 
greater weight is given in C n to those nuclear spins with 
the weakest couplings to the electron (and hence to each 
other in the effective equation of motion). Therefore 
we expect these integrals of motion to play successively 
weaker roles in the dynamics in general, and especially 
in the dynamics of the electron. Consequently, we focus 
on the C n with small positive n. 

Suppose we consider the infinite time limit for a finite 
system N < oo. Let us consider the time average of some 
physical quantity over this time, defined by 

(O) = lim - [ dtO(t) . (52) 

In ergodic systems, this limit is well-defined for most 
quantities. In integrable systems, the existence of this 
limit may depend upon which quantity is averaged. Con- 
sider the time average of Eq. (f5Tj) for n = — 1. Since 
ho = I, it is conserved and thus time-independent, and 
can be brought outside the time average. Therefore, one 
has 

(fith_i) = -I x (m) . (53) 

If we now assume both (rh) and (h_i) are well-defined, 
then the left-hand-side in Eq. (foU)) vanishes. Thus we con- 
clude I x (to) = 0, i.e. the time average of the nuclear 
field is parallel to the total angular momentum: 

(m) = «| . (54) 

The proportionality constant must satisfy \k\ < 1 since 
m is a unit vector. If |k| ^ 1, this indicates persistent 
motion of the nuclear field in the asymptotic time evolu- 
tion. Our numerics actually show that k is a convergent 
quantity. 
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D. Approximate treatment of time averages 

Despite the very large number of integrals of motion, 
a complete analytic solution of this problem has eluded 
the community. In this section, we adopt a less ambitious 
but more pragmatic approach, and attempt to calculate 
the long-time average of the nuclear field asymptotically. 
We do this by first making the quasi-ergodic assumption 
that this time average is equivalent to a uniform average 
over the accessible regions of phase space. The latter av- 
erage is analogous to a microcanonical ensemble average 
in statistical mechanics, but with more constraints pro- 
vided not only by energy conservation but by the many 
integrals of motion. Mathematically, using the integrals 
C n , one has 

(m)„ = \ fv[l 3 ]ni5(J2li-\)5(\J2 a ^\- H N) 

p j j 

p 

X [] 5 ( C n ~ Cn) , (55) 
n=2 

where c„ are the fixed values of the integrals C n , T>[Ij] = 
Y[j dlj , and we have introduced the partition function 

z p = /WiKE^-O^iE^i-H") 

3 

V 

xl[S(C n -c n ) (56) 

n=2 

In addition to this quasi-ergodic assumption, we make 
a separate approximation, which consists of truncating 
the O(N) series of constraints in Eqs.(([55] ) .([56" ]) ) to only 
the first few. The degree of approximation is deter- 
mined by the parameter p above. For p = 1, we include 
only energy and angular momentum conservation, while 
p > 1 counts the number of additional conserved quanti- 
ties taken into account. This approximation is motivated 
by the fact that, as argued above, the C n are expected to 
have increasingly little influence on the system's dynam- 
ics for increasing positive n and negative n. 

To calculate (m) p for small p, it is convenient to intro- 
duce the (unnormalized) distribution of the vector m: 

P p (m;l,H w )= fvfcWY, 1 *- 1 ) ( 5? ) 
J j 
p 

xS(^2ajIj - Hjvm) JJ 5(C n - c„) . 

j n=2 

By rotational invariance, this distribution is actually a 
function only of the angle 9 between m and I and not the 
absolute directions of either vector. Thus we may write 



with m • I = cos 9 = u. The desired expectation value is 
then obtained by a simple spherical average as 



IP ~ r i 



(59) 



J^duPpiu) 

where we have suppressed the dependence on \,hi- 

The distribution P p (cos#) can be calculated for small 
p by a saddle point technique (for N 3> 1). Details are 
given in Appendix [Bj For p = 1, one obtains the very 
simple distribution 



i\(cos(0)) = e QCOS W , 



with 



a 



3a Hjvl 



a 2 _ ~ a 2 NP ' 

and we define the moments of exchange couplings 

1 N 

a* = Ly a n 

N ^ 3 

3 = 1 

The resulting average direction from Eq.(|59p is then 

(m)i = I ( cotha 

For a -C 1, (ni)i has an even simpler expression: 

w.-lf 

The first improvement upon the above result is to in- 
clude one additional conserved quantity, which gives (see 
Appendix [Bj 



(60) 
(61) 

(62) 

(63) 
(64) 



P 2 (cos6») = e 



7 cos 0-\-(3 cos 9 



(65) 



where 



6IHMq2 • a 3 ~ aa 4 ) + js^^ ~ (« 2 ) 2 ) 
2NH N P[(a 2 Y - 2~a 2 a~a^ + (^3)2 - (J? - (a) 2 K] 
3l 2 (a^3 - (yf) 2 /[I 2 (y - a 2 )] 



=t7; (66) 



2N[(a 2 ) 3 - 2a 2 aa 3 + (a 3 ) 2 - (a 2 - (a) 2 )a 4 ] ' 

with K = 2h 2 • hi = -C 2 + 2P Y,j 9j The resulting 
average of the nuclear field direction becomes 



<iii) 2 = I 



2e 4 <3 +/3 sinh(7) 



F 



( 7+2/3 \ 
) 



F 



2v^ J 



J_ 

2/3 



P p (m;l,Hjv) = P p (cos0;l,H 



(58) 



where F{z) is the imaginary error function defined by 
F(0) = and F'{z) = -j^e z ' '. 

In principle, we can take into account more and more 
conserved quantities. However, the calculation becomes 
more and more cumbersome. 

In the next section, we compare the predictions of this 
approximate calculation with numerical simulation of the 
classical equations of motion. 
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IV. NUMERICAL SIMULATION 

In this section we briefly review the numerical simula- 
tion of the classical equations of motion for the nuclear 
spins alone. The equation we are using is Eq. flTj) . We 
use a finite difference discretization of the equations of 
motion to advance the nuclear spins from the initial con- 
ditions to any subsequent time. 

In order to protect the unit magnitude of the nuclear 
spins in the finite difference equations, we use a rotation 
matrix about the effective field Hat to advance the nu- 
clear spin states from {I,(i)} to {Ij(t + At)}. Rather 
than simply advance to the next time step, we use the 
Euler method - we take the values of li(t + At), calculate 
|(Hjy(i) + H jv(i + At)), and use this new averaged effec- 
tive field value to produce new values for Ij(t + Ai). This 
procedure is repeated a number of times until I; (t + At) 
is converged. Mathematically, this above process can be 
expressed compactly as follows, start with 

{Ut)} and H° N (t)=J2 a Mt)i 

i 

then repeat 

H k N (t + At) = J2adHt + &t), 

i 

H h N +1 (t) = ±(H° N (t) + H k N (t + At)), 

until 

|I^ +1 (£ + At) + At)| <e, 

where, the operation lZ(m, (f) is the rotation about di- 
rection m by angle 4>, the integer variable k is the times 
repeated, and e is used to measure the convergence. 

To verify correctness, we checked for conservation of 
energy and total angular momentum I for all the data 
series we generated. In our simulations we found that 
these are conserved at an accuracy of AE/E, | Al|/I ~ 
0.1% — 0.2%. Numerical simulation for a number of spe- 
cial solvable cases were also carried out ( "solvable" means 
they can be solved analytically after the short-time av- 
erage for electron spin introduced in Sec. IIII A[) . These 
consist of the problem of two nuclear spins with differ- 
ent coupling constants and the problem of many nuclear 
spins with identical coupling constants. In both cases, we 
find agreement between the simulations and the analytic 
results within the numerical accuracy. 

We compare using the Euler method with a rota- 
tion matrix to the standard finite difference 4th order 
Runge Kutta in our simulation by finding the distribu- 
tion P(cos 9) (introduced in Sec lIIIDi and measured by 
finding the histogram of the variable cos 9) . The two sim- 
ulations were performed over 500, 000 timesteps of 0.05, 
starting from the same initial conditions for 128 nuclear 
spins with coupling constants uniformly distributed from 



3500 
3000 
2500 
§ 2000 
U 1500 
1000 
500 





0.2 0.4 0.6 0.8 1 

cos(60 

FIG. 1: Comparison between histogram profiles generated by 
4th order Runge Kutta and our modified Euler method. The 
dark gray profile is from data by 4th order Runge Kutta, 
the light gray one is from our modified Euler method. As 
expected, they almost coincide with each other. 



to 1, The two measured distributions agree well, as 
plotted in FigO but the 4th order Runge Kutta method 
requires much finer time discretization (smaller At) than 
our modified Euler method in order to preserve energy 
and total angular momentum to the same degree of ac- 
curacy. Furthermore, our modified Euler method more 
naturally conserves the magnitude of each individual spin 
vector. 

In the simulation we find the complete time evolution 
of the variable cos 9. We also find the value of k defined 
in Eq.([54|). The simulations are carried out for 600,000 
timesteps of At — 0.1, with different initial conditions 
and three different coupling profiles. The approximate 
theoretical k values from Eq. (|6"3")) and Eq. (|rJ7|) are plotted 
against the numerical n values in Fig |2l3l4l for the differ- 
ent coupling profiles. For each figure, the correlations 
between numerical data and two theoretical predictions 
are explicitly calculated and listed in Table. |TJ For two 
data sets X and Y , the correlation px.Y between them 
is defined as 

cov(X,Y) 

Px.y = , (69) 

a x Oy 

where cov(X, Y) is the covariance between X and Y and 
ax,(Ty are the standard deviation of data set X and Y , 
respectively. The correlation is 1 in the case of an increas- 
ing linear relationship, —1 in the case of a decreasing lin- 
ear relationship, and some value in between in all other 
cases, indicating the degree of linear dependence between 
the variables. The closer the coefficient is to either — 1 
or 1, the stronger the correlation between the variables. 
Although for some (very few) initial configurations, the 
analytical prediction given by Eq. ([6T|) deviates from the 
numerical values more than Eq. (|63[) . the correlations in 
the table clearly indicate that analytical prediction by 
Eq.flSTD is much better than Eq.®. In addition, Eq.(lB"5D 
always gives positive k, Eq. (|67p can give negative k as 
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K data set numerics and Eg, (163[l numerics and Eg. (l67|l 

Fig|2] 0.817407 0.891817 

Fig|3] 0.771773 0.926312 

Fig[4] 0.704288 0.991488 



TABLE I: Correlation between n value data from numerics 
and two theoretical predictions given by Eq. (|63|l and Eq. (|67[) 
for three different coupling profiles, the data are the same as 
the ones used in the corresponding figures of the left column 
above. 




-0.5-0.25 0.25 0.5 0.75 1 

^numerics 

FIG. 2: Comparison of theory and numerics for values of 
K. The couplings are given by gaussian coupling profile 

dij = 2.0 exp [— 4u ^ J ' where i,j are integer values con- 
strained by i 2 +j 2 < 900/-7T, 900 nuclear spins total. The coor- 
dinate of each box(gray) is (^-numerics, fi>Eo. (|63[) ): while the 
star(black) is (ft-numerics, K-Eq. (|67|l 'l. same numerical value 
represents same initial condition. The error bar of the numer- 
ical k is too small to be seen in this figure. In this figure, 39 
different initial configurations are generated. The solid line 
is used as reference line, meaning the exact agreement be- 
tween theory and numerics. (The same conventions hold for 
Fig. EU 




numerics 



FIG. 3: Comparison of theory and numerics for values of k. 
The couplings are given by gaussian coupling profile Oi,j = 

2.0 exp [— ^/ 7r ^4oo"' ~^]; where i, j are integer values constrained 
by i 2 +j 2 < 900/tt, 900 nuclear spins total. In this figure, 30 
different initial configurations are generated. 
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FIG. 4: Comparison of theory and numerics for values of k. 
The couplings are given by gaussian coupling profile Oi,j = 

2.0 exp [— ^Yooo"^]' wnere hi are integer values constrained 
by i 2 +j 2 < 900/tt, 900 nuclear spins total. In this figure, 29 
different initial configurations are generated. 



shown by several points in Fig. I2I4| when the averaged 
effective nuclear magnetic field is anti-parallel with the 
total nuclear angular momentum I. The correlation ta- 
ble also indicates that, (though we need more numerics to 
confirm this indication) Eq. (f63|) works better for the cou- 
pling profiles which are relatively more peaked, Eq. (|63[) 
and Eq. (|67p works better for the coupling profiles which 
are more extended. Our explantion for this feature is 
that, the more peaked the coupling profile is, the better 
the approximation given by Eq. (J4j) is; and under this ap- 
proximation, the collective variables h n defined by P5|) 
can be considered roughly proportional to H^v, so Eq. (|63]) 
is probably enough to give relatively good prediction, and 
Eq. (|67[) won't improve the prediction significantly as for 
the more extended coupling profiles. 



V. ASYMPTOTIC BEHAVIOR OF (S z {t)} 

As discussed in the Introduction, even at long times 
t there will be some very weakly coupled spins which 
have not evolved appreciably. These spins, being very 
weakly coupled to the electron, also do not significantly 
contribute to its average or to the nuclear hyperfinc field. 
Thus we need consider only the dynamics of those spins 
with couplings atf > k with some constant k of 0(1). 
We will keep k as a free parameter in the considerations 
below to assess the sensitivity of our results to this choice. 
Quantities which are sensitive to O(l) changes of k are 
clearly not reliably obtained from our approximations. 

Let us for the time being focus on wavefunctions with 
the (stretched) exponential form of Eq. © , for which the 
size of the region of dynamically coupled spins evolves 
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logarithmically, as given in Eq.([7]). For such logarith- 
mic growth, the number of dynamical spins increases ex- 
tremely slowly at long times, with very few additional 
spins being added with each proportional increase (e.g. 
doubling) in time. Thus we believe that the most of 
the spins within this region have to a very good ap- 
proximation explored their accessible phase space - sub- 
ject of course to the constraints of integrability. This is 
the physical argument justifying the use of Eg. ([55)1 with 
the time-dependence only appearing through the slow in- 
crease of the set of dynamical couplings. 



In Sec. IIII C[ we argued that invariants C n have de- 
creasing relevance to the electron dynamics with increas- 
ing positive n. In the more specialized context of this 
section, we can slightly refine this contention. Specif- 
ically, consider the average (rh p ) defined in Eq.([55]) in 
the long-time limit. We are interested only in the time- 
dependence of this quantity, i.e. how it changes as addi- 
tional very weakly coupled spins lj are introduced into 
the average. These appear in the integration measure, in- 
side the definition of m, and in the <5-function constraints 
due to the constants of the motion. However, in all these 
places except the measure and the ^-function constrain- 
ing the total angular momentum I - and specifically in 
all higher invariants C n with n > 2 - the additional in- 
troduced spins are weighted by positive powers of their 
very small couplings cij. Thus we expect that all these 
approximations give very similar long-time dependence 
for p > 1. Indeed, by explicit computation (using e.g. 
Eq. (|65|) and its generalizations) we found identical re- 
sults for the leading asymptotic long-time behavior for 
p=l,2. 

We therefore discuss the simplest but sufficient case 
p = 1, in which only energy and total angular momen- 
tum conservation are taken into account. To use Eq. (|63[) 
we require a from Eq. (|61[) , which depends upon the first 
two moments of the hyperfine couplings. Consider the 
case of a Gaussian wavefunction, 7 = 2 in Eq.®. In- 
cluding those nuclear spins at radii r < R(t), one ob- 
tains a = n: rrm and a 2 — ,„°,° jta 1 where Qq gives the 



ln(t/fc) a ^ " — 21n(t/fc)' 

overall scale of the hyperfine coupling. A typical (root 
mean square ensemble average) value of the combination 

2 ~' , and for this reason the parame- 



Hjyl „ 

NI 2 

ter a 



V / ln(t/fe) : 



< 1. 



(70) 



Assuming S z (0) = 1/2, then cos6*o = ^ and we can 
evaluate the appropriate ensemble average derived in 
Eq. ([55)) . replacing the classical spin by its short-time av- 



erage and using Eg . (|39|) . (f64 

C(t) = {S z (t)) = 



1 hB 



< N (t) ' Z (^ 
2H N (t) \(t) 3 

1 

2^ 



- «2 NP 



1 



a 



6 a 2 - 

1 1 

6mT 



(71) 



where we have replaced H^/Hat with the effective 
H^(t)/Hjv(f) found within R(t) in the first line of the 
above equations, as they are very close in the asymptotic 
time limit, and in the fourth line, the ensemble averag- 
ing is taken, where [H z N (t)\ z (t)} = aNI 2 /3 is used. The 
asymptotic logarithmic decay agrees with the heuristic 
result of the introduction, and with the numerical find- 
ings of Ref. [I0IIT3I . Interestingly, the choice of k does not 
affect the prefactor of the leading logarithmic behavior 
in Eq. ([7T|) . suggesting that not only the logarithmic de- 
pendence but its prefactor might be properly obtained 
by this approach. The prefactor in (|7Tj) is indeed of the 
same order of magnitude as found in the l/S z -lnf plot 
by Ref. [l(| whose treatment is equivalent to a classical 
one with I = S = y3/2 in the equations of motion rather 
than I = S = 1/2. 

We have carried out a similar calculation for an expo- 
nential wavefunction, 7 = 1 in Eq.®. In this case, we 
obtain 



C(t) 



1 



3 m 2 (!)' 



(72) 



which is different from the asymptotic decay for the pre- 
vious example, but agrees again with the arguments in 
the introduction. This disagrees with the suggestion in 
Reffiol that the decay in this case would have the same 
form as for the Gaussian wavefunction above. However, 
no justification was given in Refill for this claim apart 
from a slow growth of l/(S z (t)} on a plot against \nt. 
Given the spread in their numerical data, it is likely that 
the form of Eq. ([72|) would produce a good fit. 

For the case of an infinitely high potential barrier in 
the quantum dot, the coupling constants vanishes at the 
boundary, and the electron only couples with a finite 
number of nuclear spin. The nuclear spins near the 
boundary have a coupling that vanishes quadratically 
with their distance from the boundary g(r) ~ (1 — r//) 2 , 
I being the radius of the confining barrier. Repeating the 
calculation above for this case, we find 



C(t) ~ ci + c 2 



(73) 



where ci and c-i are constants. This result can also be 
understood from the argument given in the Introduction. 
When R(t) reaches the boundary of the confining poten- 
tial, R(t) = l(l-y/kjt) andTV(f) = N(l- ^/kjt) 2 , where 
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N is the total number of nuclear spins. Plugging N(t) 
into Eq.© leads to the above result. 



VI. DISCUSSION 

In this article we have derived the appropriate classi- 
cal limit for the quantum central spin problem of a single 
electron spin coupled to many nuclear spins, and shown 
that this limit is indeed valid when the number of nuclei is 
large. This validates the starting point of Ref. [HI, which 
simply assumed this. Furthermore, we have made some 
approximate analytic predictions for the asymptotic be- 
havior of the central (electron) spin in this system, and 
have shown by numerical simulation that these predic- 
tions are roughly correct. Moreover, the approximations 
can be systematically improved, as we also described, 
though the technicalities involved in this improvement 
appear formidable. In the long-time limit for a finite 
number of nuclear spins, we argued that the asymptotic 
polarization of the electron spin will on average point 
in the direction of the total nuclear angular momentum 
I, regardless of the coupling constants aj. Finally, we 
have analytically derived the functional form of the de- 
cay of the electron spin polarization toward its asymp- 
totic value. This form depends on the coupling constant 
profile, and is in qualitative and semi-quantitative agree- 
ment with observations from prior numerical investiga- 
tions. Our results elucidates the physical origin of the 
long-time remanent spin polarization, and indicates that, 
remarkably, it is not connected to the integrability of the 
central spin problem. 

The quite simple explanation of the long-time spin dy- 
namics suggests that numerous extensions of these results 
should be possible, which are interesting in light of the 
relevance of the problem to experiments in quantum dots. 
Most obviously, one could likely explain the long-time 
behavior in situatiuons with substantial nuclear polar- 
ization and/or applied magnetic fields, such as treated in 
some special cases in Refs0I3. Our derivation of the clas- 
sical limit in SeclITl was general enough to apply to these 
cases, so a treatment of this problem may proceed already 
from this point. Since integrability appears not to play 
a significant role, there also appears to be no substantial 
obstacle to analysis of more complex (and non-integrable) 
situations such as two electrons in near proximity coupled 
to overlapping nuclear populations, and crossover effects 
induced by weak nuclear dipole interactions. It would 
also obviously be interesting to treat more complicated 
dynamical quantities than the remanent electron spin de- 
scribed here, such as higher order spin autocorrelations, 
electron- nuclear entanglement, and the influence of time- 
dependent fields. 

The more important and probably more formidable 
challenge is to connect the intriguing dynamical behav- 
ior of such electron-nuclear interaction problems to ex- 
periment. As was mentioned in the Introduction, the 
hyperfine interaction is the dominant source of dephas- 



ing in current single-spin experiments in GaAs quantum 
dots. A recent review can be found in Ref0. Exist- 
ing experiments use optica l 15 ' 16 or electrica l 17 ' 18 means 
to measure this dephasing. These experiments, however, 
focus upon the short-time behavior, on the order of the 
electron precession frequency, in which the nuclear field 
can be regarded as approximately static. The non-trivial 
behavior which is the focus of this paper appears only 
on longer timescales. It is possible to use spin-echo tech- 
niques to remove the trivial dephasing effect of this fast 
electron precession about the (uncontrolled) quasi-static 
nuclear fiel d 19 ' 20 . This technique (and possibly others) 
should in principle allow the unimpeded observation of 
the slower nuclear dynamics. Very recently, optical mea- 
surements have achieved the sensitivity necessary to mea- 
sure a single electron spin in a single quantum dotj 21 ' 22 
Direct observation of the electron spin in real time can 
of course also detect both short and long time dynamics. 
Given the sometimes frantic pace of experimental devel- 
opments in this field, we hope that experimentalists as 
well as theorists will make the effort to explore some of 
this very interesting physics before rushing to more com- 
plex structures. It is conceivable that the confirmation of 
our physical picture of electron-nuclear spin interactions 
could even lead to useful ideas for future devices. For 
instance, the observation that the coupling profile deter- 
mines the form of the electron spin polarization decay 
could enable control of time-dependent spin polarization 
by quantum dot shape engineering. 

VII. ACKNOWLEDGMENT 

We acknowledge discussions with V. Dobrovitski, and 
thank L. Glazman and D. Loss for stimulating our in- 
terest in the problem. This work was supported by the 
Packard Foundation and the National Science Founda- 
tion through grant DMR04-57440. 

APPENDIX A: SELF- AVERAGING OF 
ELECTRON SPIN 

In this appendix, we show that the fluctuations of the 
dynamics in initial conditions of the nuclear spin bath 
are negligible, thus demonstrating that the system is self- 
averaging. We will do so by proving Eq. (|13p . 

We first consider h = in Eq.JTU]), in which case 

<{"u>l<il^*(t)ll>K"^»<{"i}|<il^(t)l|>]{"i» . (Ai) 

where Z = (21+ 1) N , and 

>.)>.n»';!>:i ~ z 2 

(A2) 
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With Eqs. (|A1[) , (|A2p . we can express the variance of 

2 

the spin expectation value, (S z (t)) 2 — (S z (t)) as 

i E (MKll^WlDdl^WlDK^}) (A3) 



In Eq.jXll, we have {\\S z {t)\\) 2 < ± resulting in a 



variance of 



(s z (t)) 2 -{s z (t)f < ±s <MIW> = jz ( A4 ) 



so the variance is of order ~ Z = (21 + 1) . Since 

-2 



(S*(t)) ~ 0(1), we find the h -> limit of Eq.(Tn 

For the ft, 7^ case, an analogous derivation may 
be carried out, now including the Boltzmann-like factor 
H = —hj^i If as defined by Eg. (Til)]) . The variance now 
becomes 



(S»(t))* - (S'(t)) 

4 E \({rn,},\\S z {t)\\,{n t })¥e h ^ +h ^ 



{mi,rii } 



z 2 



J] ({m l }| e ^(I|5 2 (t)|i) e -«(i|^W|i)|{m. t }> • 



(A5) 



Using (i|S' z (£)|i) < I we find a bound on the variance 



(S z (t)) 2 -(S z (t)) 
1 



<4^ £<Mie _2W iM> 



1 /sinh((27+ l)ft) 



4Z 2 V sinh(ft) 

1 / tanh(|) 
4 ( v tanh((/+ \)h) 



N 



(A6) 



JV 



With (S z (t)) ~ 0(1), we have Eq.©). 

Since ~ ^"/^rfil.^ < 1, in the large TV limit, the above 

tSll.il 1 -p ) fl J 

expression decays to zero exponentially. 



significant conserved quantities. 
Z x = J dm. J PlI^^Hjv-^aJO^'-E 11 ) ' ( B1 ) 

i i 

where rh is the directional vector of Hn- First we use 
Lagrange multipliers to get the form 

Z 1 =JdraJ dxdye 4x Hjv+iy i y V [I,-] e^Ei^-^i+y- 1 *) 

(B2) 

Next we integrate over the lj variables 

Z 1 = [drhf dxd ye - H "+^'n Sina|fltX + yl) , 
J J ' J -. J - 7|ajX + y| 

(B3) 

where / is the magnitude of the nuclear spin vector. Now 
we use the thermodynamically large number of nuclear 
spins to first approximate the sum over nuclear spins by 
a distribution function p(a) for the coupling constants, 
and then to justify a saddle point approximation 

TT sin(/|aiX + y|) = &N j rfap (a) in( 3i " J ( | f J°^ l) ) 

i I \ a ^ + y\ (B4) 

~ e -^ fdap(a)(a X +y) 2 = & - Hf- px 2 +2ax-y+y 2 ) _ 

Finally we integrate over the Lagrange multipliers x, y, 
and get 

3 (Hf v +^l 2 -2aH Jv . l ) 

Z x = / drhe 2 «/ 2 <a 2 -^) . (B5) 



We define the probability distribution as follows 



Zx = / dcos(0)Pi(cos0) , (B6) 



where is the angle between rh and I. The distribution 
function then takes the form 

Pi (cos 0) = e QCOse , (B7) 

where a = fgtt^ - as defined in Sec. [TTTD] 

NI 2 (a 2 -a 2 ) 

Similarly, simply introducing another Lagrange multi- 
plier, we can find 



APPENDIX B: APPROXIMATE DISTRIBUTION 
FUNCTION FOR THE EFFECTIVE NUCLEAR 
FIELD DIRECTION 

In this appendix we find the distribution function for 
calculating (|55|) . taking into account only the two most 



P 2 (cos0) = e 7co S e+/3cos 2 e ^ 



(B8) 



where j3 and 7 are defined in Sec. IIIIDI 

For P p (cos0) when p > 2, the calculation becomes 
more and more complicated, and hard to get an ana- 
lytical expression. 
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